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Abstract. Wc prove existence, uniqueness and regularity of weak solutions of 
a coupled parabolic-elliptic model in two dimensions; we consider the standard 
equations of magnctohydrodynamics with the advective terms removed from 
the velocity equation. Despite the apparent simplicity of the model, the proof 
requires results that are at the limit of what is available, including elliptic 
regularity in L 1 and a strengthened form of the Ladyzhenskaya inequality 

<c||/||^||v/||^ 2 2 , 

which wc derive using the theory of interpolation. The model has applications 
to the method of magnetic relaxation, introduced by Moffatt (J. Fluid. Mech. 
159, 359-378, 1985), to construct stationary Eulcr flows with non-trivial topol- 
ogy- 



1. Introduction 

In this paper we prove global existence and uniqueness of solutions to the follow- 
ing coupled parabolic-elliptic system of equations related to magnctohydrodynamics 
(MHD), for a velocity field u, a magnetic field B and a pressure field p defined on 

C M 2 , as follows: 

-i>Au + Vp* = (B ■ V)B, (1.1a) 

f) Fi 

— — h (u ■ V)B - rjAB — (B ■ V)u, (1.1b) 
at 

V u = V ■ B = 0, (1.1c) 

where p* = p+ i \B\ 2 is the total pressure. Here v > is the coefficient of viscosity, 
and rj > is the coefficient of magnetic resistivity. 

This model has interesting analogies with the vorticity formulation of the 3D 
Navier-Stokes and Euler equations, as well as with the 2D surface quasigcostrophic 
equations. Recall that the vorticity formulation of the Navier-Stokes equations in 
three dimensions is 

— + (u ■ V)u> - rjAuj = (u ■ V)tt, (1.2) 

where u = K * uj is given by the Biot-Savart law, with K a homogeneous kernel of 
degree —2 (or rather of degree 1 — n in dimension n). Our two-dimensional model 
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has a very similar form — compare (1.2) with (1.1b) — but u is instead given by 

u = K * (B ® B), 

where K involves derivatives of the fundamental solution of the Stokes equation, 
and is homogeneous of degree —1. Unlike the 3D Navier-Stokes equations, for 
which existence and uniqueness of solutions for all time remains open, our two- 
dimensional model retains the essential features of the nonlinearities but admits a 
unique solution for all time. 

Indeed, the main purpose of this paper is to prove the following theorem. 

Theorem 1.1. Let ft be one of the following: 

• Q C M 2 is a Lipschitz bounded domain with Dirichlet boundary conditions; 

• n = M 2 ; or 

• = [0, l] 2 with periodic boundary conditions. 

Given an initial condition Bq G L 2 (Q) with V • Bq = 0, for any T > there exists 
a unique weak solution (it, B) of (1.1); that is, a pair of functions (u, B) such that 

u G L°°(0,T; L 2 >°°(n)) D L 2 (0,T; H 1 ^)) 

and 

B e L°°(0,T;L 2 (n)) n L 2 (0,T; H 1 ^)) 

satisfying (1.1) as an equality in L 2 (0, T; _ff -1 (f2)). Furthermore, for any T > e > 
and any fceff, 

u,B e L°°(e,T;H k (n))r\L 2 (e,T;H k+1 (n)). 

Our interest in system (1.1) arises from its connection with the method of mag- 
netic relaxation, an idea due to Moffatt [31]. He considers the related full MHD 
system: 

— + (u- V)w - vAu + Vp* = (B ■ V)B, (1.3a) 

— + (u-V)B-r]AB = (B-\7)u, (1.3b) 
V-m = V-B = 0. (1.3c) 
Formally, when r\ = 0, we obtain the standard energy estimate 

~(N|i a + ||B||i a )+f||Vu||i a =0; 

so as long as u is not identically zero, the energy should decay. Furthermore, by 
using the so-called magnetic helicity, which is preserved under the flow, we can find 
a lower bound for the energy of B: if Mm '■= J n A-B, where V x A = B is a 
vector potential for B, then 

C||B||i 2 > || J B|| 2 2 ||A|| 2 2 >^A-B^ = I^mI 2 > 0. 

In other words, the magnetic forces on a viscous non-resistive plasma should come 
to equilibrium, so that the fluid velocity u tends to zero. We are left with a steady 
magnetic field B that satisfies (B ■ V)B — Vp* = 0, which up to a change of sign 
for the pressure are the stationary Euler equations. 

These arguments are heuristic, and as yet there is no rigorous proof that this 
method will yield a stationary Euler flow. The first problem is that it has not 
yet been proved that the system (1.3) with 77 = has a unique solution for all 
time, even in two dimensions: the best known result is found in [22], where short- 
time existence of strong solutions is proved by means of a "vanishing resistivity" 
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argument, as well as a conditional regularity result; the conditional regularity result 
was later extended in [17]. 

With -q > 0, however, the existence theory for (1.3) is in a similar state to the 
Navier-Stokes equations, with global existence of weak solutions in two or three 
dimensions, and uniqueness in two dimensions; sec [15] and [36]. (Interestingly, 
global existence of weak solutions in two dimensions for the case v — but r\ > 
was proven in [23] , with various extensions in [8] and [7] , and conditional regularity 
results in [17] and [39].) 

The second problem is that, even with global existence and uniqueness, the sys- 
tem may not possess a limit state. Assuming that the equations have a smooth 
solution for all time, and furthermore that ||-B||oo < M for all time, Nunez [35] 
showed that (with r\ = 0) the kinetic energy must decay to zero, but that the mag- 
netic field may not have a weak limit when a decaying forcing / G L 2 (0, oo; L 2 (£l)) 
is added to the u equation. 

If all we are interested in is the limiting state, the dynamical model used to obtain 
that steady state is not particularly important: in a talk given at the University 
of Warwick, Moffatt [32] argued that dropping the acceleration terms from the u 
equation and working with a "Stokes" model — such as equations (1.1) — might 
prove more mathematically amenable. As a first step towards a rigorous theory 
of magnetic relaxation for this model, this paper thus establishes existence and 
uniqueness theory for (1.1) in the case r\ > in two dimensions. 

The proof of Theorem 1.1 is divided into several sections: 

• In Section 2, we introduce the weak L p spaces, denoted L p '°°, and use the 
theory of interpolation spaces to prove the following generalised version of 
the 2D Ladyzhenskaya inequality: 

II/IIl.< C ||/II^I|v/||^ 2 2 . 

• In Section 3 we consider elliptic regularity for the Stokes equations 

-vAu + Vp = V • /, 
V • u = 0, 

and show that u £ whenever / 6 L 1 . 

• In Section 4, we use the results of the previous two sections to prove global 
existence and uniqueness of weak solutions for (1.1) in a bounded domain 
Q and the whole of M 2 . 

• In Section 5 we prove higher-order estimates to show that the solutions stay 
as smooth as the initial data permits for all time, and hence that after any 
arbitrary time e > the solution is smooth. 



2. Interpolation and Ladyzhenskaya' s inequality 

In order to prove existence and uniqueness for our system (1.1), we will require 
a variant of Ladyzhenskaya's inequality. We first recall the standard inequality 
proved by Ladyzhenskaya [26]: if f2 C IR 2 is a bounded Lipschitz domain, then for 
u G H x {tt), 

HU« < cH^HI^ 2 - (2.i) 

One can prove this simply by using the embedding H 1 / 2 C L 4 and interpolating 
H 1 / 2 between L 2 and H 1 : 

NU* < 4u\\ Hl/ , < cIHI^ihi^ 2 . 
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But one can also prove it directly for u E Cl(Q) (see [26], or [18], equation (4.8) on 
p. 17): we can write u 2 = 2 J udjudxj, and then integrate u 4 = (« 2 ) 2 to get 



u 4 | dxi dx2 



< 



/max|u(a;)| 2 da;i I ( / max |w(a;)| 2 dx2 
-oo X2 ) \J-oo 11 



< 4 // maxlttc^ul dx\dx^ ) // max|u9iu| dx\dx2 



and the result follows by applying the Cauchy-Schwarz inequality The result for 
u E H x {ti) follows by density of C*(fi) in H 1 ^). 

The variant of Ladyzhenskaya's inequality that we require is the standard in- 
equality with ||u||l2 replaced with ||m||^2,oo, where L 2,00 (R 2 ) is the weak L 2 space: 

\\fh* <4f\\% 2 „ II V/H^ 2 . (2.2) 
In fact, in Section 2.2 we will prove the stronger inequality 

II/IIlp <c||/||^||/||^ p (2-3) 

for every / £ L q -°°(R n ) n BMO(R"), using the theory of interpolation spaces. 

The inequality (2.3) is not altogether new: an alternative proof is sketched in 
[25] , and it is a strengthening of the inequality 

ll/lk < cll/ll^l/H^, (2.4) 

which has been proven a number of times before; see [10], [24], [14] and [3]. In 
particular, the elegant proof of (2.4) in Chen & Zhu [10], which uses the John- 
Nircnberg inequality for functions in BMO, is adapted in McCormick ct al. [29] to 
give an alternative proof of (2.3). 

2.1. Weak L p spaces. The weak L p spaces are defined as follows (see. [20], §1.1, 
for example). Let C R™ be measurable. Given a measurable, a. e. -finite function 
/: — > M, we define its distribution function df : [0, oo) — > [0, oo] by 

df{pt) := fi{x E fl : |/(a;)| > a}. 

Then, given 1 < p < oo, the weak LP space, denoted L p '°°(f2), consists of all 
measurable, a. e. -finite functions / for which the quantity 

f C p 
\\f\\p,oo ■= inf < C > : d f (a) < — for all a > 

is finite (sec [20], Definition 1.1.5). Note that || • \\ Pl00 is not a norm, but only a quasi- 
norm — the triangle inequality fails to hold, but instead we have the replacement 
inequality 

||/+<?|k*.~ < 2(||/|| W .co + || fl ||£p.c.). 

It is fairly simple to see that L p (fi) C L p ' 0O (n): note that 

\\f\\%= I \f{x)\ P > I \f(x)\ p >aVd f (a), 

JR" J 

{x:\f(x)\>a} 

so if / E L' p (n) then d f (a) < \\f\\ p LP a~P, and hence / E L p '°°(n). However, the 
function f(x) = \x\- n / p is in LP'°°(R n ), even though it is clearly not in L p (W l ): 
notice that 



d f (a) = fj,{x E R" : |x|-"/ p > a} = fi{x E W : \x\ < a~ p/n } 

LP' 00 



where uj„ is the volume of the unit ball in R n , so 11 f 11 - — — / ' 1 ' p 
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Just as we can interpolate between L p spaces using Holder's inequality, we can 
interpolate between weak L p spaces to get to a strong LP space in between. The 
following result is proved in [20], Proposition 1.1.14. 

Lemma 2.1 (Weak-strong interpolation). Let 1 < q < p < r < oo. Then 

\\fh,<c\\f\\i-%\\f\r L ^, 

where ^ + s = I. 

q r p 

We can also use weak LP spaces to weaken the standard Young's inequality for 
convolutions; see [20], Theorem 1.2.13. 

Theorem 2.2 (Young's inequality). Let 1 < p < oo and 1 < q, r < oo satisfy 
- + 1 = - + -. Then there exists a constant C p . q . r > such that, for all f £ L p 
and all g £ L r '°° , we have 

\\f*g\\ L l-~ < C p . q . r i|/jMM|L— • 

2.2. Interpolation spaces. In order to prove our weak version of Ladyzhcnskaya's 
inequality, we will use some of the standard theory of interpolation spaces. We 
recall here the basic facts we require: for full details, see the books of Bennett and 
Sharpley [5], §5.1, and Bergh and Lofstrom [6], §3.1. 

Let (Xo, Xi) be a compatible couple of Banach spaces (that is, there is a Haus- 
dorff topological vector space X such that Xq and X\ embed continuously into X) . 
The K -functional is defined for each / £ Xo + X\ and t > by 

K(f,t) = K(f,t;X , := inf{||/o||x + : / = fo + fx] 

where the infimum is taken over all representations / = /o + /i of / with fo £ Xq 
and fx £ X\. 

For < 8 < 1, we define the interpolation space (Xq,Xi)$ )00 as the space of all 
/ £ Xo + X\ for which the functional 

||/|| e ,oo:= sup t- e K(f,t) 

0<t<oo 

is finite. A very useful property of interpolation spaces is the estimate on the norms: 

ll/lkoo^cll/lCll/lliv (2.5) 

when < 6 < 1 (see [6], §3.5, p. 49). 

As a simple example of interpolation, note that 

(L 1 (W l ) , L°° (R n )) 1 _ 1 / p oo = L p '°°(R n ) 

if p > 1 (see [5], Chapter 5, Theorem 1.9). This is hardly surprising given that 
the definition of interpolation spaces is modelled on that of weak L p spaces. In 
fact, this equality remains true with L°° replaced with BMO (see [5], Chapter 5, 
Theorem 8.11): 

(L 1 (R n ), BM0(R")) 1 _ 1 / PiOO = L p '°°(R n ). 

The so-called reiteration theorem allows us to interpolate between interpolation 
spaces: it says that when we interpolate between two interpolation spaces of the 
same couple (Xq, Xi), we get another interpolation space in the same family. 

Theorem 2.3 (Reiteration Theorem). Let (Xq,Xi) be a compatible couple of 
Banach spaces, and let < 9q < 0\ < 1- Set Aq = (Ao,Ai)e o 00 and A\ = 
(Jr ,^i)ei,oo. If0<e< 1, then 

{Ao, Ai)e,oo = (Xo, Xi)o' :00 

providing & = (1 - 6)9 Q + 66\. 
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The proof may be found in [5], Chapter 5, Theorem 2.4, or [6], Theorem 3.5.3. 
Using this, we can prove a weak version of our generalised Ladyzhenskaya inequality 
(2.3). 

Lemma 2.4 (Weak interpolation). For any f G L q '°°(K n ) n BMO(K"), and any 
q < p < oo, 

\\fh^<c\\f\\iu\f\\B^ p - 

Proof. By [5], Chapter 5, Theorem 8.11, we have 

L*°°(R n ) = (L 1 (R"),BM0(]R n )) 1 _ 1/g , oo 

provided that 1 < q < oo. Set 33 := (L 1 (R 7l ),BM0(R n ))i )0O) and note that by 
(2.5) we have ||/||<b < C||/||bmo- By the Reiteration Theorem (Theorem 2.3), we 
obtain 

L p '°°(R n ) = (L q '°°(M. n ), 03) QiOO 
with q < p < oo, provided that a solves 1 — i = (1 — a)(l — i) + a ■ 1, or in other 
words that a = 1 — q/p. Thus, using (2.5), we obtain 

as required. □ 

By combining this with Lemma 2.1, we obtain our generalised Ladyzhenskaya 
inequality (2.3). 

Lemma 2.5 (Strong interpolation). For any f G L 9 '°°(R") n BMO(K"), and any 
q < p < oo, 

||/lk.<c||/||i£. 

Proof. Given p > q, choose any r and s such that q<r<p<s< oo. Then 

||/||£-<c||/||^||/||£.,oc, 

where a satisfies + ^ = i. Applying Lemma 2.4 to each of the two factors on 
the right, we obtain 

Wfh* < cCcii/iii^ii/ii^^CcII/iii^ii/iIbm^) 

<r M\ f \\( 1 - a )i/ r + a ^/ s \\ fii( 1 - Q )( 1 -'?/ I ')+ Q ( 1 -9/ s ) 
-i HI J \\li<<*> \\j IIbmo 

= 011/11^-11/11^', 
as required. □ 

It follows from Poincare's inequality that W 1,n C BMO in n > 2 dimensions 
(see §5.8.1 in [16]): so, in particular, in two dimensions H 1 C BMO. (One may 
also prove that H n ' 2 c BMO in n dimensions: see Theorem 1.48 in [4].) Thus for 
/ e L 2 '°°(R 2 ) n ij^M 2 ), setting n = 2, p = 4 and q = 2 in Lemma 2.5 we obtain 
(2.2): 

\\fh< <c\\f\\%*„\\Vf\\% 2 . 

When is a bounded Lipschitz domain in K 2 , we may extend a function / G Hq (0) 
by zero outside fl and apply the above inequality on R 2 to obtain the same for il. 

When = [0,1] 2 with periodic boundary conditions, however, a different ar- 
gument using Fourier series is required: we obtain the same inequality using the 
Sobolev embedding L 4 C H 1 ^ 2 and the fact that 

/ = E f^ ik ' x => ii/Hl*<ck 1/2 ii/iil*.~ (2.6) 

|fe|<K 
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(a weak form of Bernstein's inequality; see [33, 34] for the relevant theory of Fourier 
series and McCormick et al. [29] for the proof). Indeed, writing 

J = / y Ike + 2^ fo e 

\k\<K \k\>K 

we obtain, using (2.6) and L 4 C H 1 / 2 , 

||/IU<<CK 1/2 ||/|U> 




<CK 1/2 \\fh 
<CK 1/2 \\f\\l 

Minimising over n we obtain (2.2). A similar argument involving Fourier transforms 
can be used to obtain a more general version of (2.2) and (2.3) on the whole space; 
see McCormick et al. [29]. 

3. The Stokes operator and elliptic regularity in L 1 

We now consider the Stokes equation alone. Take /: fl — > R 2x2 , and define 
V ■ / : ft — > M 2 componentwise as follows: 

2 
i=l 

With abuse of notation, we write L v {yi) = L p (£l;M. 2x2 ) where no confusion can 
arise. We consider the equations 

-vAu + Vp = V • /, (3.1a) 
V-« = 0, (3.1b) 

with Dirichlet boundary conditions if 51 is bounded. By setting / = B g) B 
(i-C fi,j = BiBj) we recover equation (1.1a): 

2 2 

(V • (B ® B))j = di(BiBj) = BAB, = [(B ■ V)B] 3 , 

i=l i=l 

as B is divergence-free. In this case, if B € L 2 (f2), then B <E> B is in L 1 (J7), so 
the right-hand side behaves like the derivative of an L 1 function. If / e L p (Vl) for 
p > 1, one would expect that u S W 1,p (£l), but this docs not hold for p = 1. If 
it did, in two dimensions we would obtain u € VF 1,1 (f2) C L 2 (fl). In fact, in this 
section we prove that, when / G L 1 (r2) in (3.1), then u € L 2 '°°(il). 

The solution of this equation is given by convolution with the Green's function 
or fundamental solution of the equation: let U ', g solve 

-vMJ + Vq = 5, 

w-u = o, 

where 8 denotes the Dirac delta function. Then the solution of (3.1) is given by 

u = U*(V-f), p = q*(V-f). 

Integrating by parts, if duU € L 2 '°°(£l), then by Young's inequality (Theorem 2.2) 
we can show that 1 1 1 1 j^a, < C||/||^i. 
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In the case fl = R 2 , we have explicit formulae for the fundamental solution U 
and q: 

1 



U ^ X) Am, 



S - S v log \x\ 



/ \ ^ Xj 

(A derivation of the fundamental solution may be found in §IV.2 of [19].) Integrat- 
ing by parts with respect to k, we obtain 

Ui (x) = [U * (V ■ fMx) 

2 2 



V / d k U itj (x-y)fk,j(y)dy. 



Now, 



and so 



d k Ui ! j(x) = — 



&ikXj ~\~ &kjXi XiXjXfc ^ Xfc 



\dkUi,j(x)\ < 



TTU\X\ 

As noted in Section 2.1, ^ is in L 2 '°°(R 2 ), and |[<9 fc t/ M || L 2,oo < Using Young's 

inequality (Theorem 2.2), we obtain 

||u|U a .« < c\\d k U id \\ L ,,oo\\f\\ L , < c\\f\\ L i. (3.2) 

Thus, whenever / e L^M 2 ), u € L 2 '°°(R 2 ). 

In the case where £1 = [0, l] 2 with periodic boundary conditions, one can also 
write down an explicit formula for the fundamental solution — see [21] and [11], for 
example — and obtain (3.2) again; the details are very similar to the above case, 
and we omit them. 

In the case where f2 is a bounded Lipschitz domain, while we no longer have 
an explicit formula for the Green's function C7, by Theorem 7.1 in [30] we have 
V(7 G L 2 oo (n) whenever f! C R 2 is a bounded Lipschitz domain, thus using 
Young's inequality (Theorem 2.2), we obtain (3.2) on a bounded Lipschitz domain 
as well; i.e., whenever / £ i 1 (^), u e L 2 '°°(Q). 

4. Existence and uniqueness of weak solutions 
We return now to the system 

-i/Au + Vp* = (B ■ V)JB, (4.1a) 

— + (u ■ V)B - f/AB = (B ■ V)m, (4.1b) 
V-m = V'B = 0, (4.1c) 

where = p + \ |-B| 2 - We will show that equations (4.1) have a unique weak 
solution for all time in the three cases of f2 described in Theorem 1.1. We first 
define a weak solution. 

Definition 4.1. A weak solution of (4.1) on [Q,T] is a pair of functions (u,B) 
such that 

u e L°°(0,r;L 2 '°°(n)) n^fO.T;^ 1 ^)) 
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and 

B e L oo (0,T;L 2 {Q))nL 2 (0,T;H 1 (n)) 
satisfying (4.1) as an equality in L 2 (0, T; H" 1 ^)). 

Note that the pressure p is uniquely determined by u and B by solution of a 
standard elliptic boundary value problem; see [9] or [18]. We will prove the following 
theorem. 

Theorem 4.2. Given B 6 L 2 (Q) with V • B = 0, for any T > there exists a 
unique weak solution (u,B) of (4.1). 

In Section 4.1, we will prove existence of a weak solution in the case fi C R 2 is a 
Lipschitz bounded domain with Dirichlet boundary conditions, while in Section 4.2 
we prove existence of a weak solution in the case H = R 2 . The proof of existence 
in the case where Q = [0, l] 2 with periodic boundary conditions is analogous to the 
previous two, and we omit it. Finally, in Section 4.3, we prove uniqueness of weak 
solutions. 

4.1. Global existence of solutions in a bounded domain. In this subsection 
we prove existence of a weak solution on a Lipschitz bounded domain f2 C R 2 , with 
Dirichlet boundary conditions, using the method of Galerkin approximations. To 
do so, we first set up some notation. Let H := {u 6 L 2 (S1) : V • u = 0}, and let 
II be the Leray projection II: L 2 (il) — > H, i.e. the orthogonal projection from L 2 
onto H . 

We define the Stokes operator as A := — IIA. Let {4> m }meN C C°°(f2) be the 
collection of eigenfunctions of the Stokes operator on fl with Dirichlet boundary 
conditions, ordered such that the eigenvalues associated to (f> m are non-decreasing 
with respect to to. Let V m be the subspace of H spanned by <j)\, . . . , <fi m , and let 
P m : H — > V m be the orthogonal projection on to V m . 

In order to use the Galerkin method, we consider the equations 

-vAu m + VpT = {B m ■ V)B"\ (4.2a) 

— + P m [{u m ■ V)£? m ] - r,AB m = P m [(B m ■ V)w m ] , (4.2b) 
V • u m = V • B m = 0. (4.2c) 

Thinking of u m as a function of B m given by equation (4.2a), it is easy to check that 
(4.2b) is a locally Lipschitz ODE on the finite-dimensional space V m , and thus by 
existence and uniqueness theory for finite-dimensional ODEs (Picard's theorem), 
there exists a unique solution B m e V m of equation (4.2b), with u m given by 
equation (4.2a). 

Proposition 4.3 (Energy estimates). The Galerkin approximations are uniformly 
bounded in the following senses: 

u m g L^iO^T; L 2 '°°(n)) D L 2 {0,T; H 1 ^)); 

B m g L°°(0,T;L 2 (n)) n L 2 (0, T; If 1 (f2)). 

Proof. Take the inner product of equation (4.2a) with u m and the inner product 
of equation (4.2b) with B m , and add to obtain 

~\\B m (t)\\ 2 L2 +^1^^(^11^+771^5-^11^=0. 
Integrating over [0, t] we obtain 

\\B m (t)\\ 2 L2 +2v f \\Vu m (s)\\ 2 L2 ds + 2r l [ \\ VB m (s) ||| 2 ds 
Jo Jo 

= \\B m (0)\\ 2 L2 <\\B \\h, 
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and taking the essential supremum over all t G [0,T], it follows that 

esssup||B m (t)||| 2 +2u f \\ V« m (s)|| 2 a ds + 2n f || VB m (s)||| 2 ds 
te[o,T] Jo Jo 

<2||Bo|li 2 , 

which proves that 

u m G L 2 (0, T; tf^Q)); B m G L°°(0, T; L 2 (0)) n £ 2 (0, T; if^Q)). 

As in Section 3, the solution u m to equation (4.2a) is given by convolution with 
U, the Green's function for the Stokes equations. By (3.2), we have 

||u m (t)IU a .~ <c\\(B m {t)f\\ L , <c||B m (t)||i a , 

so taking the essential supremum over t G (0, T) tells us that 

u rn g ^(O.TjL 2 ' 00 ^)), 

which completes the proof. □ 

Proposition 4.4. 77ie Galerkin approximations are uniformly bounded as follows: 

OB 11 



G L J (0 } T;H- L (SI)). 



dt 

Proof. Taking the H^ 1 norm of the B equation yields 
dB r 



dt 



< V \\B m \\ H i + \\P m [(B m ■ V)u m }\\ H -> + \\P m [(u m ■ V)B m ]\\ H - 1 . 

H- 1 

To estimate the nonlinear terms, we take the inner product with a test function 
4> G -ffg(^) with = 1 and then take the supremum over such </>: 

\\P m [(B m ■ V)« m ]|| = sup \{P m [(B m ■ V)u m ],cf>}\ 

= sup \{{B m -V)P m ^u m )\ 

<SMV\\B m \\ L ,\\u m \\ L 4V^\\ L , 

< \\B™\\ L 4u m \\ L *. 

Similarly \\P m [(u m ■ V)B m ] < ||B m || L 4 ||M m || L 4. By applying Ladyzhenskaya's 
inequality (2.1) to B m and our weak Ladyzhenskaya's inequality (2.2) to M m , we 
obtain the the following estimate: 



dB m 



<r;|| J B™||^+ C ||B m !| i .||B m || Hl ||«™|| i ^||« m || ffl 

H- 1 



dt 

as required. □ 
So, in summary, we have the following: 

u m are uniformly bounded in L°°(0, T; L 2 '°°(f2)) t~l £ 2 (0, T; i/ 1 (Sl)), 
B m arc uniformly bounded in L°°(0, T; L 2 (fl)) n L 2 (0, T; _ff 1 (f2)), 

arc uniformly bounded in L 2 (0, T; H^ 1 ^)). 

To extract a convergent subsequence of B m , we use the Aubin-Lions compact- 
ness lemma: 

Theorem 4.5 (Aubin-Lions compactness lemma). Let X C B C Y be Banach 
spaces such that the inclusion X C B is a compact embedding. Then, for any 
1 < p < oo and any 1 < q < oo, the space 

f : / G LP(0,T;X) and % G L^(Q,T;Y) 
at 
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is compactly embedded in L p (0,T; B). 

Proof. The original result of Aubin [2] and Lions [27] covers the case when 1 < p,q < 
oo. Chapter 3 of Temam [38] contains both the original case (see Theorem 2.1, 
p. 185), as well as the case p = 2, q = 1, whenever X and Y are Hilbert spaces (see 
Theorem 2.3, p. 187). The general case (and many other similar results) is proved 
in the paper of Simon [37], §8, Theorem 5 and Corollary 4. □ 

We use Theorem 4.5 together with the Banach-Alaoglu compactness theorem to 
extract a subsequence, which we relabel as B m , such that 

B m ^B in L°°(0,T;L 2 (ft)), 

B m — B in L 2 (0,T;ir 1 (fi)) ) 

B m -> B in L 2 (0,T;L 2 (O)), 

T^f ^<o,r;tf-<n)). 

Since the limit B G L°°(0, T; L 2 (£7))nL 2 (0, T; /^(fi)), it is straightforward to show 
that (B ■ V)B G L 2 (0,T;If _1 (n)). This allows us to de/ine u to be the unique 
solution of 

-vAu + Vp* = (B ■ V)B, (4.3a) 
V ■ u = 0, (4.3b) 

where u G L°°(0, T; J L 2 '°°(fJ)) n L 2 (Q, T; ii" 1 ^)) by standard elliptic theory for the 
Stokes equations (see Section 3 above, and Theorem 2.3 in Chapter 1 of [38]). 
Having defined such a u, we now want to show that u m does indeed converge to 
this u in the appropriate senses; this will allow us to show that the nonlinear terms 
involving u converge and thus that the B equation is satisfied in the limit. 

Proposition 4.6. u m — > u strongly in L 2 (0, T; L 2 '°°(J7)), and u m — * u weakly in 
L 2 {0,T;H l {Q)). 

Proof. Subtracting the equations for it m and u, we obtain 

-vA(u m - u) + V(pJT - P*) = V • (B m 8> B m -B ®B) 

= V • [J5 m <g> (J5 m — B) + (B m -B)®B]. 

By elliptic regularity, we obtain 

\\u m - u\\ L 2,^ < c\\B m <g> (B m - B)|| L i + c||(B m - B) ® B|| L i 

< c||B m || L3 ||S m - B|| £a + c||B™ - B\\ L2 \\B\\ L , 

= c\\B™-B\\ L z{\\B™\\ L * + \\B\\ L *) 

<c(K + M)\\B m -B\\ L *, 

where K = sup meN sup te r 0T i ||B m || i 2, M = sup te r 0T i ||B||i2. Squaring and inte- 
grating in time yields 

f \\u m {t) u(i)\\l^ dt<c[ T \\B m (t) - B{t)\\l- 2 At. 
Jo Jo 

As the right-hand side converges to zero, so is the left-hand side, and hence u m — s- it 

strongly in L 2 (0, T; L 2 '°°(f2)). Let v be the weak limit of u m in L 2 (0, T; fl^fi)); 

it remains to show that it = v. As iJ 1 (il) C L 2 (Q) C L 2 '°°(fi), we have 

(£ 2 <°°)*(f2) C L 2 (0) C H-^fl). So if « m v in L 2 (0, T; i? l (ft)), then u m r> 

in L 2 (Q, T; L 2 '°°(f2)) (because we are testing with a smaller set of functionals). 

But it" 1 — > it strongly (and hence also weakly) in L 2 (0,T; L 2 ca (fl)), and thus by 

uniqueness of weak limits u = v, and the proposition is proved. □ 



12 



D. S. MCCORMICK, J. C. ROBINSON, AND J. L. RODRIGO 



We now proceed to show that the nonlinear terms in the B equation converge. 
The following proposition is symmetric in B and it, and thus applies to both the 
(it ■ V)B and (B ■ V)it terms. 

Proposition 4.7. Suppose that: 

• u m -)• it and B m -)• B (strongly) in L 2 (0, T; L 2 <°°(£1)); and 

• u m , B m are uniformly bounded in L°°(0, T; L 2 '°°(f2)) fl L 2 (0, T; H 1 ^)). 

Then (after passing to a subsequence) 

P m [{u m -V)B m ] {u-V)B in L 2 (0, T; 

Proof. We begin by showing that P m [(it m ■ V).B m ] are uniformly bounded in 
L 2 {0,T;H-\VL)): observe that 

\\P m [(u m ■ V)B™]|| H -i = sup|<(it m ■ V)(P m 0,B TO )| 

< IKMlB^BupllWlUa 

< II^II^IIB-H^ 

<c||« m ||^||V M m ||^ 2 ||B'"!|^||VB'"||^ 2 

(where the supremum is taken over all <j> G J? 1 (fi) with = 1). Thus, squar- 

ing and integrating in time and applying Holder's inequality shows that P m [(u m ■ 
V)B TO ] is uniformly bounded in L 2 (Q, T; H" 1 ^)), and hence that a subsequence of 
P m [(u m ■ V)B m ] converges weakly-* in L 2 (0, T; PT _1 (fi)); as usual we relabel this 
subsequence as the original sequence. 

To show that the limit is indeed (it • V)J5, we test with a slightly more regular 
test function. Let <f> G C°(0, T; H 1 ^)). Then 

[ (P m [(u m ■ V)B m ] - (it • V)B, 0) dt 
Jo 

= [ (P m [(u m -S7)B m -(u- V)B],d>)dt+ f ((u-V)B,P m 0-$dt. 
Jo Jo 



For the second integral, note that 

II < \\(u-V)B\\ H -i\\P m <l>-<l>\\ H i -+0 
as m — > co. For the first integral, we have 

1= f ((it m • \7)B m - (it • V)J3, P m <j>) dt 
Jo 

= f <(it m ■ \7){B m -B) + ((it m - it) • V)B, P m <f>) dt 
Jo 

< / (||lt" l || i 4||B" l -B|| i 4 + ||lt" l -lt|| i 4||B|| i 4)||V0|| L2 dt 

JO 

< max ||V0|| i2 / T (||tt m ||i4||B ro -B|U4 + ||tt ro -u|| L 4||B|| i 0d*. 
te[o,T] Jo 
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Now, by the weak Ladyzhenskaya inequality (2.2), 
T \\u m \\ L 4B m - B\\ L , dt 







< / ||u'"t' 2 .l|Vi,"'||^||iJ™-iJ||^.l|V(B"'-B)t/ 2 

^(firt,.) I/4 .(f iiv«-iii.) 1/4 - 



|B m -B||i a 



o 



^ T ||V(B--B)||^ 



The first, second and fourth integrals arc bounded independent of m, and the third 
tends to zero as m — > oo, so f Q T ||it m || i 4 \\B m — B\\ L 4 dt — > as m — > oo. Similarly, 
J T \\u m — u\\ L 4\\B\\ L 4 dt — > as m —> oo. Hence 

/ (P m [(t6 m -V)B m ]-(M-V)B,0)dt^O for all £ C°(0, T; iJ 1 (f2)). 
Jo 

Thus P m [(u m ■ \7)B m ] ^ (u ■ V)B in L 2 (0, T; by uniqueness of weak-* 

limits. □ 

Hence (u,B) is indeed a weak solution of (4.1), which completes the proof of 
Theorem 4.2 in the case where H is a Lipschitz bounded domain in M 2 . 

4.2. Global existence of weak solutions in R 2 . We turn now to the proof of 
Theorem 4.2 in the case H = R 2 . Instead of Galerkin approximations, we mollify 
the equations, and then show convergence as e — > 0. The arguments, though, are 
not so different from those in the previous section, so we only outline the main 
changes. 

Let p > be a smooth function with compact support (i.e. p S C^°(R™)) such 
that J Rn p = 1. Given e > 0, we define p £ by p £ (x) :~ -^p{x/e). and we define the 
operator J £ by J £ v = p e * v, i.e. 

(J e v)(x) = / p s (x - y)v(y)dy. 

For brevity we will write v := J £ v. We consider the mollified MHD equations on 
the whole of R 2 as follows: 

-i/Am £ + Vpl = J £ {(& ■ V)£? £ ] , (4.4a) 
dB e _ ~ 

— + J £ [(ut ■ V)B £ ] - i]J £ ABt = J £ [(Bc • V)w e ], (4.4b) 
V ■ u e = V ■ B £ = 0. (4.4c) 

As in the previous section, we think of u e as a function of B e . Then, using standard 
properties of the mollifier J £ (see [28], Lemmas 3.5 and 3.6), it is straightforward to 
show that equation (4.4b) is a locally Lipschitz ODE on if — one basically follows 
the proof of Proposition 3.6 in [28]. Hence by the Picard theorem for infinite- 
dimensional ODEs, equation (4.4b) will have a unique solution B e € H, so long as 
||JB e || L 2 remains finite, and u £ is given by equation (4.4a). 
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Repeating the estimates of Propositions 4.3 and 4.4, with slight modifications to 
account for the extra mollifiers, we again have the following: 

u £ are uniformly bounded in L°°(0, T; L 2 '°°(R 2 )) n L 2 (0, T; ij^R 2 )), 

B E arc uniformly bounded in L°°(0, T; L 2 (R 2 )), 

J e B e arc uniformly bounded in L 2 (0, T; ij^R 2 )), 

dB e 

-gj- are uniformly bounded in L 2 (0, T; i? _1 (E 2 )). 

Because we are working on R 2 , we cannot apply the Aubin-Lions compactness 
lemma (Theorem 4.5) directly (because the embedding H 1 C L 2 is no longer 
compact). Instead, there exists a subsequence of B s which converges strongly 
in L 2 (0, T; L 2 {K)) for any compact subset K C R 2 (see Proposition 2.7 in [9]), and 
the limit satisfies 

B G L°°(0,T;L 2 (R 2 )) nL 2 (0,T;iJ 1 (R 2 )). 

Thus, we may again define u to be the unique solution of equation (4.3). A 
modification of Proposition 4.6 shows that a subsequence of u 6 converges strongly 
to u in L 2 (0, T; L 2 {K)) for any compact subset K C R 2 , and this local strong 
convergence allows us to pass to the limit in the nonlinear terms: an argument 
similar to Proposition 4.7 will show that (after passing to a subsequence) 

J e [{u* • V)B £ ] (u ■ V)B, J e [(B e • V)m £ ] a (B • V)u 

in L 2 (0, T; 7f _1 (R 2 )) (see §2.2.4 of [9] for full details). Thus (u, B) is indeed a weak 
solution of (4.1), which completes the proof of Theorem 4.2 in the case f2 = R 2 . 

4.3. Uniqueness. We now prove that weak solutions are unique. Note that the 
following proof applies equally in all three cases of Theorem 1.1. 

Proposition 4.8. Let (uj,Bj), j = 1,2, be two weak solutions with the same 
initial condition Bj(0) = Bq, such that 

Uj G L°°(0, T; L 2 >°°(n)) n L 2 (0, T; H 1 ^)), 

B 3 e L oo (0,T;L 2 (il))nL 2 (0,T;H 1 (n)). 

Then iti = 112 and B\ = B-2 as functions in the above spaces. 

Proof. Take the equations for (ui,Bi) and (112,-62) an d subtract: writing w = 
«i - «2, z = B\ — B2 and q — pi — P2, we obtain 

-vAw + \7q= (J5i • V)z + (z ■ V)J5 2 , (4.5a) 

dz 

— + (ui ■ V)z + (w ■ V)B 2 - r]Az = (B 1 ■ V)w + (z ■ V)u 2 , (4.5b) 
at 

V ■ w = V • z = 0. (4.5c) 
Taking the inner product of (4.5a) with w yields 

H|VH| 2 2 < (HBill^ + \\B 2 \\ Li )\\z\\ L 4Vw\\ L z, 
so dividing through and using Ladyzhenskaya's inequality we get 
^||Vw|U 2 < llzll^dlBill^ + HBall^) 

|^4 2 ||V*||^ 2 (||^i||^ 2 ||V J Bi||^ 2 -I- ||^ 2 ||^ 2 ||V-B 2 ||i4 2 ) 
l^llV^/aiVB!^ 2 + ||VB 2 ||^ 2 2 ), (4.6) 

since Bj G L°°(0, T; L 2 (il)). By the elliptic regularity arguments from Section 3, 
we obtain 

\\w\\ L2 ,oo < c{\\B x z\\ L , + ||B 2 z|| L i) < c||*|| £a (||Bi|| ia + \\B 2 \\ L *) < c\\z\\ L 2. (4.7) 
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Using the weak Ladyzhcnskaya inequality (2.2), we obtain bounds in L as follows: 

\\w\W < ^INII^UVzll^dlVBxIl^ 4 + IIV^H^ 4 ). (4.8) 

By a similar argument, we can take the inner product with Uj of the original 
equation for Ui (1.1a) and obtain 

\\Vu 3 \\ L 2 < cWBjWlt < cWBjWvWVBjWv, 
so using (2.2) we get 

KHl* < cll^ll^HV^II^ 2 < cll^ll^HVB,^/ 2 < (4.9) 

using elliptic regularity. 

Since §f e L 2 (0, T; and z g L 2 (0, T; ff 1 ^)), we can take the the inner 

product of (4.5b) with z and obtain 

^NlL + r/llv* 

<c||V«||L a ||tl»|| i 4(||B 1 || £ 4 + || J B 2 || L 4)+ C ||V^|| L 2||^|| L 4(||W 1 || L 4 + ||u 2 M. 

I II 

We bound I and II separately, using Ladyzhenskaya's inequality several times. By 
using (4.8) and Young's inequality with e — \ for 8/5 and 8/3 on I, we obtain 

I < l\\zf^\\Vz\\%\\\VB l \\% i + HVB^I^^dlVBill^ 2 + IIV^H^ 2 ) 

< ^\\z\\% i \\Vz\\%\\\VB 1 \\% A + ||VB 2 || 3 / 2 4 ) 

Similarly, by using (4.9) and Young's inequality with e = \ for 4/3 and 4 on II, we 
obtain 

II < l^ll^llV^ll^dlV^II^ 2 + ||VB 2 || 4/ 2 2 ) 
< \\\Wzf L2 + ^llzlH.dlVBrlli, + ||VB 2 ||| 2 ). 

Hence 

^Mh + v\\V*\\h < -^\\z\\U\\VB 1 \\l 2 + \\VB 2 \\l 2 ). (4.10) 
By Gronwall's inequality, for all t < T we have 

||*(i)||| 2 < ll^olll, exp ^||VB 1 ( S )||| 2 d S + ^ t ||VB 2 ( S )|| 2 2 d s 

As Bj G L 2 (0, T; 7? 1 (r2)), the integrals in the exponential are bounded, and since 
the initial conditions were the same, Zq = 0, and hence ||2|| L 2 = on [0,T]. By 
(4.10), ||Vz|| L 2 =0on [0,T]. Hence, by (4.7), \\w\\ L *,oc = on [0,T], and by (4.6), 
|| Vto|| £ 2 = on [0, T}. This proves uniqueness of weak solutions. □ 

This completes the proof of Theorem 4.2. 



5. Higher-order regularity estimates 

In this section, we prove the second part of Theorem 1.1; that is, that the solution 
(it, JB) becomes smooth after an arbitrarily short time e > 0. In particular, we prove 
that if we start with initial data in H k (H), then the solution stays in H k (H.) for all 
time: 
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Theorem 5.1. Let k e N. Suppose B e H k (n) with V • B = 0. Then, for any 
T > 0, the unique weak solution of (4.1) satisfies 

u,Be L°°(0,T;H k (n))nL 2 (0,T;H k+1 (n)). 

This immediately implies that the solution (it, B) becomes smooth after an ar- 
bitrarily short time e > 0. 

Corollary 5.2. Given any T > e > and any fc G N 7 the unique weak solution of 
(4.1) satisfies 

u,Be L°°(e,T;H k (n)). 

Proof. Fix e > 0. We already know that u,B G L 2 (0, T; iJ x (ft), so for some time 
h < e/2, u(<i),B(<i) G if 1 (ft). Applying Theorem 5.1, we obtain 

u,B G L°°(e/2,T;iJ 1 (ft)) n L 2 (e/2, T; iJ 2 (ft)). 

Furthermore, if we know that 

m,JB G L°°(e(l -2 1 - fc ),T;ff fc - 1 (ft)) ni 2 (e(l -2 1 - fc ),T;ff fc (ft)), 

then there is some time <fc such that e(l — 2 1_fe ) < < e(l— 2 _fc ) and u(tk), B(tk) G 
H k (fl), and so applying Theorem 5.1, we obtain 

u,B G i°°(e(l - 2- fc ),T;ff fc (ft)) ni 2 (e(l -2- fe ),T;i7 fc+1 (ft)), 

The result follows by induction on fc. □ 

We will prove Theorem 5.1 by induction on fc. Since the base case and the 
induction step require rather different arguments, we split them into two separate 
propositions. 

Proposition 5.3. Suppose B G iJ 1 (ft) with V ■ B = 0. Then, for any T > 0, 
the unique weak solution of (4.1) satisfies 

u,B G L°°(0,T; H 1 ^)) n L 2 (0,T; H 2 (n)). 

Proof. Take the inner product of (4.1a) with —Au, the inner product of (4.1b) 
with —AB, and add: 

l_d 

2dt' 

= ((u • V)J3, AB) - ((B • V)u, AS) - ((£? • V)J3, Au). 

Now 

|((u- V)B,AB>| < || U |U4||VB|U4||AB|| i2 

^cllull^UVnll^llVBll^llABll^ 2 
^^IIABIIi.+cllwlll^llV^lll.UVBII 2 , 

and 

|((B-V)«,AB)| < llBH^IIVull^llABll^ 

^cIlVtill^llAull^llflll^llVflll^llABlU, 
< ^||A£?|| 2 2 + ^||A«||| 2 +c||V«|| 2 2 ||B|| 2 2 ||VB|| 2 2 

and 

|((B-V)B,A«)| < IIBII^IIVBII^HAmIU. 

< C ||B||^ 2 ||VB|| i2 ||AB||^ 2 ||Atx|U 2 
< J||Au||i, + |||AB||i a + c||B||i a ||VB||i a 



--llVBlH.+HIAnlli.+^IIAiJII 2 . 
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hence 

~l|Vfl||i, + ^||A«||i, + 5||Afl||i, 

< C ||VB|| 2 L2 (||«||i a .-||VTi||i a + ||V«|| 2 L2 ||B||i 2 + IIBIIMIVBHIO . (5.1) 

Since the integral of the last bracket is finite, by Gronwall's inequality we get 
that B £ L°°(0, T; H 1 ^)), and hence reusing this bound in (5.1) yields u,B £ 
L 2 (0,T; H 2 (il)). Finally, take the inner product of (4.1a) with u to obtain 

v\\Vu\\l, < ||B||£«||V«|[ i3| 

so 

v\\Vu\\v < \\B\\l, <c||B|| La ||VB|| L 2, 
and since the right-hand side is bounded, u £ L°°(0, T; H 1 ^)). □ 

To prove the induction step, we will need a higher-order estimate on the nonlinear 
term. 

Lemma 5.4. Let s > n/2 be an integer, and let u £ H s (ft) and v £ H S+1 (SY) such 
that V • u = V • v = 0. Then 

\\(u-V)v\\ H s < Cv\\u\\ H s\\v\\ H s + r. 

Proof. When f2 is a bounded domain, this follows easily from the fact that H s is 
a Banach algebra for s > n/2 (see Theorem 4.39 in [1]). For the periodic and W™ 
cases, this actually holds for any real number s > n/2; see Lemma 10.4 in [12]. □ 

With this in hand, we proceed to the proof of the induction step. 

Proposition 5.5. Fix an integer k > 2. Let Bo £ H k (D,) with V-Bo = 0. Suppose 
that 

u,B £ L°°(0, T; H k ~ 1 (Cl)) n L 2 (0, T; H k (fl)). 

Then 

u,B£ L°°(0,T;H k (n)) n L 2 (0, T; B fe+1 (fi)). 

Proof. Take the inner product of (4.1a) with (— l) k A k u, the inner product of (4.1b) 
with (-l) fc A fc B, and add: 

~\\B\\lk + i^\\u\\ 2 Hk+1 + n\\B\\ 2 Hk+1 

= (-l) k [{(B • V)u, A k B) + ((B • V)B, A k u) - ((« • V)B, A fc B)] 

< c [\\(B ■ V)u\\ Hk \\B\\ h h + \\{B ■ V)B\\ m \\u\\ Hk +\\(u- V)B\\ Hk \\B\\ H „] 

< c[\\B\\ 2 Hk \\u\\ Hk+1 +2\\B\\ Hk \\B\\ Hk+1 \\u\\ Hk ] 

< + l\\B\\W + c + \\B\\ 2 Hk ) \\B\\% k 

so 

f t \\B\\ 2 Hk + v\\u\\ 2 Hk+1 + V \\B\\ 2 Hk+1 < c\\B\\ 2 Hk (\\u\\ 2 Hk + \\B\\ 2 Hk ) . (5.2) 

Since the integral of the last bracket is finite, by Gronwall's inequality we get 
that B £ L°°(0,T;H k (Q,)), and hence reusing this bound in (5.2) yields u, B £ 
L 2 (0,T;H k+1 (n)). 

Finally, if k > 3 take the inner product of (4.1a) with (— l) fc_1 A ^ to obtain 
v\\uf Hk <c||(B.V)B|| Hfc -i||u|| H *-i < llBll^HBll^llull^-i, 
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and since the right-hand side is bounded, u <E L°°(0, T; H k (fl)). In the case k = 2, 
Lemma 5.4 does not apply to \\(B ■ V)B\\ H k-i , and so instead we take the inner 
product of (4.1a) with —Aw and estimate as follows: 

v\\Auf L2 < \((B-V)B,Au)\ < ||B|U«||VBM|Au|U a , 

so 

\\Au\\ L2 KWBW^WVBhiWABWyi 
and since the right-hand side is bounded, u g L°°(0, T; H 2 (n)). □ 

This completes the proof of Theorem 5.1, and hence also Theorem 1.1. 

6. NON-RESISTIVE CASE (?] = 0) 

In the above we have developed an essentially complete theory of existence, 
uniqueness, and regularity for the system (1.1) when rj > 0. 

The non-resistive case (rj = 0) is much more difficult, and analogous to the 
vorticity formulation of the 3D Eulcr equations in the same way that the resistive 
system has similarlities with the 3D Navier-Stokes system (as discussed in the 
introduction). Two-dimensional models with similar structure to these canonical 
3D equations (such as the 2D surface quasigeostrophic equation [13]) have attracted 
considerable attention in recent years, and we plan to present an analysis of (1.1) 
with rj = in a future paper. 
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